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AFIT/GSO/ENY/02-4 


Abstract 


Previous efforts have been directed at the guidance and control of free flying 
satellites clusters using reaction thrusters. A Tethered Formation of Satellites with 
distributed sensors has great potential to enhance surveillance and imaging of earth objects. 
To maintain the shape of the formation and keep the tethers in tension, the system needs to 
be spinning. To take the advantages of distributed sensors, the spin axis of the formation 
must have a component on nadir direction. General study has been focused on a planar 
formation and the results showed that Earth pointing configurations are not stable. This 
study demonstrates that spin stabilization and the use of gravity gradient for formation 
flight with tethered satellites reduce the necessity of thrusters for station keeping 
maneuvers. A tethered satellite formation, which consists of three satellites at the corners 
of an equilateral triangle, and two end bodies at the opposite side of this triangular plane is 
studied to obtain the stability based on the Likins-Pringle equilibrium for rigid satellites, 
for long period. Parameters that affect the stability of the tethered formation and the 
application of the Likins-Pringle relative equilibrium to the tethered formation are 
presented in this study. Equilibrium conditions for the tethered formation and the stable 
formation design for a long period are demonstrated. Depending on the size of the 
formation and the mass ratio of the satellites to end bodies, long-term stability is 


achievable for a continuously earth-pointing tethered satellite formation. 
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STABILITY OF A TETHERED SATELLITE FORMATION 


ABOUT THE LIKINS-PRINGLE EQUILIBRIA 


I. Introduction 


1.1 Motivation 


Using a satellite cluster for a space mission has become an interesting subject in 
recent years. Using several small satellites in a cluster instead of a large satellite has 
numerous advantages. For some space missions such as surveillance or terrain mapping, 
in order to improve resolution, the sensors must have a large aperture, which directly 
increases the size and the mass of the satellite. On the other hand the size of the satellite 
is limited by the capability of the launch vehicles. However with distributed sensors, it 
may be possible to achieve a resolution that cannot be attained with a single satellite. 

In some applications, each small satellite in the cluster communicates with the 
others and shares the processing, information and mission functions. A cluster of this 
kind might be used to form a large aperture. While current developments may allow 
apertures the size of meters, it may be possible to form an aperture with the size of 
kilometers with distributed satellites. This feature can be useful for space based radar and 
remote sensing missions, which require large aperture products to achieve acceptable area 


coverage and precision or lower resolution [1]. 


1.1 


Satellites formations have received a lot of attention in recent years and have been 
the subject of many research efforts. The U.S. Air Force is exploring formations for use 
in surveillance, passive radiometry, terrain mapping, navigation, communications and 
ground target identification as mentioned in the research of Irvin [2]. 

Unfortunately to keep the group of satellites in relatively close proximity in a 
formation may require significant fuel expenditure [3]. Station keeping maneuvers and 
burns to keep the formation are due to the perturbations and unmodeled forces. An 
alternate concept is to form a formation with small satellites, using tethers between these 
satellites in the formation [4]. Instead of using thrusters for each satellite for all station 


keeping maneuvers, using a tether between satellites would reduce the fuel consumption. 


1.2 Problem Statement 


Although using tethers between satellites in the formation is beneficial to reduce 
the fuel consumption, a stable configuration must be found in order to make the idea 
practical. All the satellites in the formation have different forces than a single satellite in 
orbit. While the main force is the gravitational force on the satellites, there are also 
tension forces through the tether on the satellites in the formation. To keep the formation 
in a desired configuration, tethers between the satellites must be in tension. These tension 
forces will help to rigidize the tethered satellite system and keep the formation from 
collapsing or otherwise altering its shape. 

The center of mass of the formation follows an orbital path that is assumed to be 


elliptical. However satellites in the tethered formation have affects on each other with the 
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tension force through the tethers between them. Thus the dynamics of a tethered satellite 
formation are affected by the tensional forces. These additional forces on a satellite cause 
the satellite to behave in a non-Keplerian manner [5] [6]. 

Since the studied application of the tethered formation is Earth surveillance and 
sensing, the system must be designed so the sub satellites are widely distributed about the 
nadir directions. Previous studies have been focused on a planar formation and the results 
showed that Earth pointing configuration is not stable. 

For the stability of the tethered formation, to use the benefits of the gravity 
gradient forces, two end bodies are placed at the opposite side of the triangular plane of 


the satellites as in Figure 1.1. 






l.end body 


Sat 3 


Sat 2 


2.end body 


Fig 1.1 Tethered Formation 


1.3 


1.3 Objectives 


The goal of this thesis is to demonstrate a stable tethered satellite formation for a 
long period. This thesis focused on a closed-loop formation with two end bodies at the 
opposite sides of the triangular plane of the satellites as shown in Figl.1. Since the 
application of the tethered formation is Earth surveillance and sensing, the system will be 
designed so the sub satellites will be widely distributed about the nadir directions. To find 
the stability condition and the equilibrium position of the tethered formation about the 
nadir direction, analysis of Likins-Pringle relative equilibrium and the possibility of using 
this equilibrium will be addressed. 

Stability of the formation will be analyzed for a range of spin rates, formation 
sizes and formation masses. Due to the complexity of the dynamics, the equations of 
motion will be solved numerically. Parameters that affect the stability of the formation 


will be presented with the results. 


1.4 


II. Background and Previous Research 


In this chapter benefits of using tethers in space will be presented with some 
possible applications. Research about tethered satellite systems and the dynamics and 


modeling of tethers will be presented. 


2.1 Possible Space Missions With Tethers 


Over the last two decades, engineers have realized the potential benefits of tethers 
in space. It has been pointed out that it is possible to use a tether between different bodies 
for many space missions. 

Using tethers between distributed satellites to form a large aperture is a possible 
way of increasing the capability of a surveillance system. The physical size of an 
aperture, which directly relates to the achievable resolution, is limited by the launch 
vehicle capacity. On the other hand, the properties of an earth surveillance system can be 
increased with distributed small satellite connected with tethers. 

Creating artificial gravity with the rotation of two bodies connected with tethers is 
an important issue for long duration manned space flights. To separate an orbital station 
into two parts connected by a tether and spin the system about its center of mass is the 
basic idea of creating artificial gravity. As it is mentioned in reference Dynamics of 


Space Tether Systems, Chobotov was the first investigator to present a detailed dynamic 
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analysis of this motion and Gemini 11 was the first spacecraft to demonstrate this 
concept in 1966, which connected to the rocket stage Agena with tether [7]. 

Docking of the orbiter to a remote tethered unit, which is separated with a long 
tether from the station, is one possible way to avoid docking impacts on the station and 
to save considerable amount of fuel during maneuvering in this step [7]. 

An orbiting interferometer, which has different receivers connected by tethers 
could be beneficial to provide high resolution. 

Scientific measurements in the atmosphere at altitudes of about 100 km could be 
possible with tethered units, which are connected to the orbiter at higher altitude. This 
altitude of atmosphere is not accessible for aircraft because of the low density and also is 
not accessible for the spacecraft because of the high density for space missions at the 
same time. However an atmospheric probe, which is connected to the orbiter with a very 
long tether can solve this problem [7]. 

A tethered unit, which can be deployed from a spacecraft into a dense atmospheric 
layer can transfer the aerodynamic breaking force to the vehicle through the tether. After 
the desired reentry conditions are reached, the tethered unit can be separated from the 
orbiter. During this maneuver it is possible to save a considerable amount of fuel [7] [8]. 

Another useful application of tethers in space is the electrodynamic tether system, 
which has a conducting tether deployed from the station. Interaction between the tether 
and the geomagnetic field of earth could supply enough power for orbital maneuvering 


of the station [7]. 
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2.2 Previous Research 


There has been a great deal of research about tethers in space over the past two 
decades. Most of this research is related to the dynamics and modeling of the tethers and 


tethered satellite systems. This research will be presented in two different sections. 


2.2.1 Tether Modeling 

In order to design missions and predict the behaviour of the tethered satellite 
system, engineers had to develop some techniques to model the tether. While most of the 
researchers have considered massless tether for the satellite formation, some engineers 
have modeled the tether with different techniques. To model the tethers, which are used 
in tethered satellite system, some methods must be used to discritize the continuous 
tether. There are some basic differences between discrete and continuous systems. 
Primarily, discrete systems have a finite number of degrees of freedom, and continuous 
systems have an infinite number of degrees of freedom. Discrete systems are governed 
by ordinary differential equations while continuous systems by partial differential 
equations [9]. The difficulties in obtaining closed-form solutions to problems of 
continuous systems can be summarized as the difficulty of solving the partial differential 
equations of the continuous systems. So to get closed-form solutions for a tether, it 
should be discritized by some techniques. Approximate solutions for the continuous 
systems generally involve discritization of the systems, and these approximate solutions 
differ from each other with the discritization technique that is used [10]. 

The most commonly used techniques to discretize the continuous tether are the 


assumed mode method and discritization with particles (bead model). Discritization with 
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particles method models the tether with small masses connected with massless tethers. 


All the external forces can be applied to these small particles and it is possible to obtain 


ordinary differential equations. 





Figure 2.1 Bead Model for Tether [7] 


Forces on the particles can be shown as: 


Vist 





Ax. alee Ax ——> 


Figure 2.2 Forces on the Tether [7] 


The assumed modes method discretizes the continuous system prior to the 
derivations of the equations of motion. The solution for the system is assumed in a finite 


series form of space-dependent functions, but the coefficients are time-dependent 
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generalized coordinates instead of being constant. This series is substituted in a kinetic 
and potential energy, thus reducing them to discrete form, and the equations of motion 
are derived by means of Lagrange’s equations [10]. 

Biesbroek et. al.[11] has compared the results of two different modeling methods 
of tethers. These two methods are the bead model and the continuum tether model. A 
comparison of two computer simulation models has been given by the author. These two 
simulations are TETSIM (bead model) and DATES (continuum tether). In this report it is 
pointed out that when a tether is modeled with beads, which are connected by 
inextensible and massless lines, it is possible to obtain accurate results for many 
applications. And also it is possible to reduce the simulation time with bead model, even 
if many beads are used. 

Modi et. al.[12] have studied the dynamics and control of tethered spacecraft 
during deployment and retrieval. In this research, the author has examined some of the 
important aspects of the other studies about tethered systems, which include the 
modeling of tethers and system dynamics and control. 

Dynamic behavior of a tethered system in elliptical orbits has been investigated by 
Takeichi et. al.[13]. For this study, the tethered system is modeled with a combination of 
rigid bodies, and the effects of tether elasticity and atmospheric drag is studied. 

Dynamical characteristics of a tethered stabilized satellite have been researched by 
Ashenberg et.al.[14]. It is pointed out by the author that, attaching a tether with a mass at 
the opposite of the satellites would increase the stability of the satellite. In this work, 


stability characteristics of a tethered satellite have been presented. 
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Zhu et. al.[14] has focused on the dynamical model and the mathematical 
formulations of tethered systems with two satellites on a line. He presented two 
dynamical models and formulations in his report for the tethered satellite system, which 
has a subsatellite and one main satellite with a point mass much greater than the 
subsatellite. For the first model, an elastic continuum is used to model the tether with 
mass and he has assumed the tether would remain straight. In the second model tether 
between the main satellites and subsatellites is modeled with discrete masses, which are 
connected with massless springs. For the second model the main satellite has much 
greater mass than the subsatellite, so he assumed that the center of mass of the tethered 


system is located at the main satellite. 


2.2.2 Multi- Satellites 


Although most research has focused on the tethered satellite systems with two end 
bodies, there are some others which have focused on the tethered systems with 
multibodies such as Kalantzis et. al. [15]. In this report the equations of motion for a 
multi-satellite tethered system in a Keplerian orbit are derived. His model considers a 
multi-satellite system which the satellites in the formation are connected to each other by 
flexible tethers. Equations for the model have been derived using the Lagrangian 
approach, and a continuum model has been used for tethers. Assumed mode methods is 
used to discritize the deformations of the tethers. 

DeCou[16] has focused on an orbiting stellar interferometer which has three 
distributed satellites at the corners of an equilateral triangle. This satellite formation 


system is determined that it is spinning around the axis, which is perpendicular to the 
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formation plane. Centrifugal force created by the spinning of the formation keeps the 
tether taut. For this system, gravity gradient, solar radiation pressure and thermal 
expansion forces are taken into account as a perturbation forces. Tethers between the 
satellites are determined under the influence of the centrifugal force and modeled with 
small point masses. Solution for the static problem has been obtained by deriving the 
nonlinear differential equations, which relate the position of each tether point to the 
tension at each point. Numerical solutions for these equations are presented. For this 
study it is assumed that three satellites at the corners of the equilateral triangle follow a 
circular path by spinning the formation around the axis perpendicular to the plane. The 
spin rate for DeCou’s formation is much greater than the orbital mean motion, so the spin 
axis maintains a nearly inertially fixed direction. 

Tragesser [4] has focused on a satellite ring as DeCou. Tethers are used to reduce 
propulsive station keeping maneuvers to maintain the formation. Earth surveillance is the 
desired application for this project, so the tethered satellite system is designed so the 
relative spin axis of the formation has a component along nadir direction. Unlike 
DeCou’s work, Tragesser’s formation spins at a rate comparable to the orbital motion. 
Conical equilibrium has been chosen to get the maximum projected area onto the earth 
for surveillance besides stability. Satellite formation has composed of three satellites at 
the corners of an equilateral triangle and the tethers between these satellites have been 
modeled with an arbitrary number of point masses. A common law of elasticity has been 
used to find the tension force on the tether. Likins-Pringle relative equilibria [17] for the 


conical case has been used to find the angle between the spin axis and orbit normal to 
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maximize the projected area onto the earth. As a result of this paper it is observed that 
the rigid body equilibrium doesn’t give the same stability results for flexible systems. 

Williams and Moore have examined the feasibility of a multi-tether system that is 
very similar to the one presented here. They focused on tethered satellite formations 
which have the spin axis aligned with the nadir direction. To use the gravity gradient 
effect for stability, they have studied a configuration with four satellites on the plane 
perpendicular to the spin axis and two anchor satellites at the opposite sides of this plane. 
This tethered satellite configuration also has a central tether between these two anchor 
satellites. In this research all the tethers are assumed as massless and remain straight 
during the motion. Clohessy-Wiltshire equations are used to describe the equations of 
motion of the satellites with respect to the center of mass of the entire system. The 
authors presented a stable configuration with small periodic oscillations of two degrees 
about the velocity vector of the reference orbit [20]. 

This thesis builds on this previous work to show that the nadir-oriented formation 
of Williams and Moore is actually an approximation of the Likins-Pringle relative 
equilibria for rigid bodies. Using this theory, this thesis demonstrates that there exists a 


range of stable configurations within a couple degrees of nadir. 
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Ill. Methodology 


In this chapter, modeling of the tethered satellite system will be presented. 
Orientation of the TSS will be given and for the stability of the formation system, spin 
stabilization on orbit will be explained. Likins-Pringle relative equilibrium concept for 
the spinning rigid body will be revisited and it will be adapted to the tethered satellite 


system. 


3.1 Nominal Design 


The formation design for this thesis is comprised of three satellites at the corners of an 
equilateral triangle and two end masses at the opposite side of this triangular plane with 
the same distance from the triangular plane. The three satellites and the two end bodies 
are considered as point masses. Tethers are placed between the satellites and between the 
end bodies and satellites with the same properties as shown in Figure 3.1 in the body 
fixed reference frame. 

Because the formation is axisymmetrical, the center of mass of the system is at 
the center of the triangular plane. In the Figure 3.1, L represents the distance between two 
end bodies and D is used for the distance between the satellites. All the tethers in the 
formation are assumed as extensible and remain straight during the motion. Modeling of 


the tethers will be discussed in Chapter 4 with more detail. 


Sal 


1. end body 





2. end body 


Figure 3.1 TSS Model 


The three satellites are at the corners of an equilateral triangle can be shown in 


Figure 3.2. 
Sat 2 
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we 


Sat 3 
; Sat 1 
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Figure 3.2 Satellites on triangular plane 
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where c is the distance between the satellites and the center of mass and b is equals to 4 


of the distance between satellites. 
Positions of the satellites and the end bodies with respect to the center of mass in 

the body frame F, can be written by geometry of the system as: 

(a) *=[0 -1/2 0 

(b) 75=[0 L/2 0] 

(c) F=LD/2 0 DiQy3)| 

(d) rslo 0 -p/V3| 

(e) F=|D/2 0 DA2V3)| (3.1) 
where 7; and 73 denotes the two end bodies, and 7’, 7.3 and 7.° denotes three satellites. 


sl? “52 53 


By using these equations, the distance of each mass to the center of mass can be found as: 
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= 1/2 (3.2) 
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7.3) = |F3| = D/v3 (3.3) 














For the orientation of the formation in body frame the inertias of each mass with 


respect to the center of mass are: 


L/4 0 

I,=m,,| 0 0 0 (3.4) 
0 O L/4 
V/4 0 

Lm) 0" 90> 6 (3.5) 
0 O L’/4 


a2 





is =m, 0 
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0 (3.6) 





(3.7) 


0 (3.8) 
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where m,, and m,, denotes the mass of the end bodies and m,, denotes the mass of the 


satellites. For the TSS model mass of the satellites are equal to each other, and so the 


mass of the satellites can be denoted as m, and the mass of the end bodies can be denoted 


as m,. Therefore the total inertia of the system is: 
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0 
0 
mL’ 
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D?m. 





I 0 0 
=|0 I, 0 (3.9) 
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The tethered satellite system, which consists of distributed satellite on a 


equilateral triangle, must rotate about the axis perpendicular to the plane to keep the 


tether in tension and to rigidize the system. Without spinning it is not possible to maintain 
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the stability of the satellite formation with flexible tethers. It is possible to rigidize the 
system by spinning the formation about the axis perpendicular to the plane of the 


satellites. 


There are three types of equilibria for the spinning rigid body, these are 
cylindrical, hyperbolic and conical equilibria (Fig 3.3). For the cylindrical and hyperbolic 
cases, because of the direction of the spin axis, the satellites’ path forms a line when 
projected onto the Earth’s surface. These formations do not take advantage of the 
distribution of sensors. On the other hand for the conical case, the spin axis has a 
component in the nadir direction, so the projection of the formation enhances Earth 


sensing capability. 





(a) Thomson Equilibrium - Cylindrical Case —(b) Likins- Pringle Equilibrium Hyperbolic Case —_(c) Likins-Pringle Equilibrium ~ Conical Case 


Figure 3.3 Types of relative equilibrium for a spinning symmetrical satellite [17] 


This thesis will focus on the conical type of equilibria to increase the projection of 
the satellite on the Earth’s surface. In the following sections Likins-Pringle relative 
equilibrium concept will be represented and then adapted to the satellite formation with 


flexible tethers. 
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For an axisymmetrical satellite with axial and transverse inertias 7, and J, 


respectively, the simplest case is to maintain the spin axis nominally normal to the orbit, 


without any gravity gradient torque on the satellite [17]. 


3.1.1 Thomson Equilibrium 
An explanation of the Thomson Equilibrium will help to explain the Likins- 
Pringle conical equilibrium, used in the formation design. For this equilibrium, the 


satellite is assumed an axisymmetrical rigid body with the transverse and axial inertias J, 


and J, respectively. The reference motion for the Thomson equilibrium is in Figure 3.4 


NE 


pin 


C 





Figure 3.4 Thomson equilibrium [17] 


For this reference motion four frames are used, which are the body frame, orbital 


frame, inertial frame, and reference spin frame. 
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ran : Body fixed principal-axis frame 


Ff. : Orbital frame 
FF : Inertial frame 
F : Reference spin frame 


The reference spin frame, F’,, rotates about the a, axis with respect to F, at rate 


v . Transformation between the reference spin frame and inertial frame is: 
cos(v‘) O —sin(Vt) 
RR? SRW)S|\ 0 1 0 (3.10) 
sin(vt) O  cos(Vt) 
where R“ indicates the transformation between the reference spin frame F’,and F, , and 
R, indicates the rotation matrix about / th axis. 


For the reference motion, body frame Ff, and the reference spin frame F, are 


aligned, but with small excursions @. The linearized transformation between these 


frames is: 


R™=|-a@, 1 a (3.11) 
3 1 


The relationship between the orbital frame F’, and the inertial frame F; is: 
cos(a@.t) 0 sin(@.t) 
R° =R,(-@,t) = 0 1 0 (3.12) 


—sin(@,t) 0 cos(@.t) 


where @, is the nominal motion. 
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By using these relationships, transformation between the body frame F, and the 
orbital frame F’, can be calculated: 
R’ = RY ROR’ SRR (RO) (3.13) 
From equation (3.10) and (3.12) 
R“(R°)" = R,(Vt)R,(@,t) = R,[(V +, Jt] (3.14) 
Substitution (3.14) into (3.13) yields: 


1 a, —-a,||cos(Vv+@,)t) 0 —sin(v+@,)t) 
R” =|-a, 1 a, 0 1 0 (3.15) 
CO. 2G 1 sin(V+@,)t) 0 cos((V+@,)t) 


The reference satellite is an axisymmetrical rigid body. The axial inertia is 7, and 
transverse inertia is /,, and because of the symmetry /,=/,=J/, and J, =J,. The 


angular velocity of this satellite with respect to the inertial frame can be written using the 
relations between the frames. The body frame is aligned with the reference spin frame 


with small excursions and the reference frame rotates about the a, axis with respect to 


the inertial frame. 


oO" =o" +0" (3.16) 
For small @, this can be approximated by: 
@” = a,b, +b, + a,b, +Va, (3.17) 


The last term of this equation can be transformed into the body frame by using 


the transformation between body frame and reference spin frame as: 
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0 1 a, -a,||0 
va,=R™\v\|=|-a, 1 @ |\v (3.18) 
0 a. -@ 1 0 
Va, = 0t,vb, + vb, —a,vb, (3.19) 


By substituting the equation (3.19) into equation (3.17), the angular velocity of the 
satellite in body frame is: 


bi 


o” =o! b,+a% b, +0" b, (3.20) 


@" = Gb, + a,b, + a,b, +01,vb, + vb, — a,vb, (3.21) 
where @” indicates the components of the angular velocity for b, axis. Equating 
components of (3.20) and (3.21) gives: 

(i) 0," =a, +a 

(ii) o,” =, +V 

(iii) @," =a,-aVv (3.22) 

To find the equations of motion it is needed to calculate the gravity gradient 

torque for the reference satellites. Gravity gradient torque depends on the position of the 
satellite with respect to the inertial frame and the inertias of the satellite. Gravity gradient 


torque of the satellite in the body frame is: 
M=M,b,+M,b,+M,jb, (3.23) 
where M, is a scalar value and indicates the magnitude of the gravity gradient torque 


for the b, axis. The components of the gravity gradient torque are given by [18]: 
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(i) M, =4yzu,-1,) 
r 
- 3ul 
(ii) M, ==4xzU1,-1,) 
r 


(ii) M, = *#yycr, =A) (3.24) 
r 
where X, Y and Z are the position of the satellite in body frame with respect to the 
inertial frame: 
7’ =Xb, + Yb, + Zb, (3.25) 


These coordinates can be found with the transformation between the orbital frame 


and the body frame as: 
hho aR 0 (3.26) 


1 a, -—a,||\cos(Vv+@,)t) 0 —sin(v+@,)t) || 0 
r°=|-a, 1 a, 0 1 0 0 (3.27) 
a, -a, 1 sin(v+@.)t) 0 cos(Vv+@,)t) || r 


—rsin(V + @,)t—@, cos(V + @. )t 
r° =| arsin(v+@,)t+@rcos(v +a, )t 


= 


(3.28) 


jo) > >> 


Ww 


—Q,rsin(v+@,)t+rcos(v +a, )t 
Therefore the coordinates X, Y and Z are: 

(i) X=-rsin(v+@,)t—@, cos(V+@, )t 

(i) Y=a,rsin(v+@,)t+a@,rcos(V+@, )t 


(iii) Z=-a,rsin(v+a,)t+rcos(v +a, )t (3.29) 


Substituting these values into the gravity gradient equations and ignoring the 


higher order terms yields: 
M, = s r? (I, —1, 1a, sin[(v + @, )t]cos[(v + @, )t] + @, cos?[(v + @, )]} (3.30) 
M, =0 (3.31) 
M, = “t =r? (I, —1, ea, sin*[V + @, )f]— a, sin[(v + @, )t]cos[(v + @, )t]} (3.32) 


For a circular reference orbit the first term can be simplified as: 


ov =5=0, (3.33) 
The equations of motion can be written [18]: 
(1) 1,0," '+(,- -1,)o;'o% =M, 
(ii) 1,0," +, -1, ook =M, 
(iii) 1,@%+(U, -1,o"o; =M, (3.34) 


where the term @” indicates the derivation of the angular velocity and can be written by 
taking the derivatives of the equations (3.22). 
(i) 0,” =&,+a,Vv 
Gi) 0," =a, 
(iii) @,” =&,-ayv (3.35) 
Substituting the gravity gradient torques and the derivatives of the angular 


velocity into equations (3.34), (Appendix B) the linearized equations of motion can be 


written as: 


&,+(1-k, var, +k,v’a, +30@°k, (a,c + a,s)c =0 (3.36) 





di, =0 (3.37) 
a, —(1-k, va, +k,v’Q, +3@7k, (a,c + a,5)s =0 (3.38) 
where: 
s=sin(v+@.)t c=cos(V+@. )t oO. = ra and 
I,-I 
k, = a (3.39) 


The linearized equations of motion include time variable coefficients because of 
the nominal motion of the satellite and the reference spin rate. It is possible to simplify 
these equations by transforming them to a new set of coordinates that will eliminate the 
periodic coefficients. The new parameters are defined by Hughes as[17]: 


(1) 4,0) =a@,(t)cos(v +a, )t+a,(t)sin(v + @, )t 


(ii) 7, (t)=a, (0) 


(il) 7;()=-—a@,(t)sin(v+ @,)t+a@,(t)cos(v + @, )t (3.40) 


Substituting these new parameters into the equations (3.36) through (3.38), the 


equations of motion can be written as: 


¥,-(V+20, +k,v)y, +[k, 80; -Vo,)-@,(V+@,)]7, = 0 (3.41) 
7, =0 (3.42) 
¥,+(Vt+20,+k,v)y, -—[k,vo,+0.(V+@,)|7; =9 (3.43) 


These three equations now have constant coefficients. The new variables y,, 7,, 
y,can be explained as three small angles between the orbital plane F, and a new frame 


F., which is assumed identical with orbital frame for the reference motion. 


Transformation between the orbital frame and the new frame, which is called 


stroboscopic frame by Hughes, is [17]: 


1 3 —% 
R“=|-y, 1 ¥, (3.44) 
Yr eae ai 1 


Because all perturbations occur between these two frames, transformation 
between the body fixed frame and stroboscopic frame is: 
R”® =R[(vV +a, )t] (3.45) 
It is easy to see that R” R” must be equal to R™R": 
RR = RR? (3.46) 
The transformation between the stroboscopic frame and the inertial frame is: 
Roa R (3.47) 
And the transformation between reference spin frame and the inertial frame is: 
R“ =R,(t) (3.48) 
Substituting the equations (3.47) and (3.48) into the equation (3.46) gives: 
R,[(V+@, )t]R* R,(-@,)= R™ R, (Vt) (3.49) 
Transformation between stroboscopic frame and the orbital frame becomes: 


R“ =R,[-(v+@, )t]R™ R,[V +o, yt] (3.50) 


As a result the transformation matrix between stroboscopic frame and the orbital 


frame can be written as: 


1 Vo VG c 0 sjj 1 Go =O, |e: Oe 8 
= «il y, |j=| 0 1 Oj;/-a@, 1 a, 0 1 O 
Y, = ol = 5 0. se) @: =a: s 0 ¢ 


where: 


c=cos[((V+q@,)t] and  s=sin[(v+q.t)] 


Multiplying out equation (3.51), gives: 


1 Vo ZY, c+s? a,c-A,s —a,(c? +s’) 
-y, 1 Y, |=|—-a,c+a,s +s? A,8+Q,c 
Vs. =H A a(c?+s°) -a,st+a,c +s? 


where: 
cos*[(v+@, )t]+sin’[(v + @, )t]=1 
Equating elements of each matrix on both side of equation (3.53) gives: 
(i) 4,0) =a@,()cos(v+@a, )t+a,(t)sinv +a, )t 
(11 ) 7, (t)=@,(t) 
(iii) 4¥;()=-a@,(d)sin(v+@,)t+a@,(t)cos(v+ a, )t 


which agrees with equations (3.40). 


(3.51) 


(3.52) 


(3.53) 


(3.54) 


(3.55) 


3.1.2 The Likins-Pringle Relative Equilibria 
The Likins-Pringle relative equilibria extend the Thomson equilibria to the case 
where the spin axis is not normal to the orbital plane [17]. The term “relative 
equilibrium” denotes a spin axis that remains fixed with respect to the orbital frame F.. 
Likins-Pringle relative equilibria allow a non zero torque unlike the Thomson equilibria, 
where the spin axis is normal to the orbit and the gravitational torque disappears. If the 


stroboscopic frame F’, is defined as the frame intermediate to the body frame F, and the 


orbital frame F'., the transformation between these frames is: 
R™ =R”R” and R” =R[(V+a,)t] (3.56) 
In the Thomson equilibrium for the reference motion without any perturbation, 
R* is equal to unit matrix. For the Likins-Pringle relative equilibrium it is a constant 
matrix, but not equal to the unit matrix. 


RY =R* R(®) (3.57) 
The term R,* indicates the small excursion angles between the stroboscopic frame 


and the orbital frame, and ® is the rotation angle about the forward direction of the 
orbital frame ¢,. 
The absolute angular velocity of the satellite is: 


wo” - o” a oo” (3.58) 
wo” =(v+a,)b, + RO" (3.59) 


The angular velocity of the satellite in stroboscopic frame is: 


Oo” =R,[-(v+o,)]O” +0" (3.60) 


The absolute angular velocity of the stroboscopic frame F’ is: 
Oe (3.61) 
Equation (3.61) includes the nominal motion and the small excursion angles. The 
term @* indicates the angular velocity due to the excursion angles. For small excursion 
angles 7, equation (3.61) can be approximated by: 
DO" =7,5, +8, +48; - OC, (3.62) 
The last term of the equation (3.62) can be transformed into the stroboscopic 
frame by using the transformation between the orbital frame F and the stroboscopic 


frame F’ as: 


—@,C,=R*| -—O, |=- @, c; (3.63) 


c Cc 


where the term c; is the second column of the rotation matrix R”. 
For the equilibrium condition without any perturbation: 


o* =0 (3.64) 


Ill 


Therefore the absolute angular velocity of stroboscopic frame is: 
@" =— QO, Ciy 8, —@, Cy $4 —, Cy 8 (3.65) 
where the terms c;, indicates the elements of the rotation matrix R*. 
Substituting equation (3.65) into the equation (3.60), the angular velocity of the 


satellite in stroboscopic frame becomes: 


DO” =O, §,+@,8,+0,8, (3.66) 


where the scalar terms in the equation (3.66): 
(i) 0, =—O, Cy 
(ii) @,=-@,c;,+(V+@O.) 
( iii ) 0, =- 0, C3 (3.67) 
To formulate the equations of motion in the stroboscopic frame F,, gravity 
gradient torque in this frame can be found by the equation for the symmetrical 
satellite[ 17]: 


M,=M,5,+M,8,+M,5, (3.68) 


U, = LT )e53€33 
M =30? 0 (3.69) 
CL, —1,)e1303 


Equations of motion become: 
(i) 10," +(U,-1,)o%o% =M, 
(ii) I,@," +(, -I, oo? =M, 
(iii) 1% +(U, -1,)o"o} =M, (3.70) 


Differentiation of the angular velocity can be calculated with the equation: 


Siew Pee : 
0" =—O" +O" x0" (3.71) 


dt 


Thus the derivation of the angular velocity of the satellite in stroboscopic frame is: 


0" =0,C3,(V+,)8,—O, C,(V+O,)5; (3.72) 


substituting equations (3.67), (3.69) and (3.72) into equations (3.70), the nontrivial 


components of the equation of motion are: 
(i) k,(3c53¢33 —¢5,€3,) + +k, )U+vV)c3, =0 
(1) k:6e..c, = 6 JF Orkut Ye, =0 (3.73) 
where k,=(/,—-J,)/I, and V=v/q@,. These equations are satisfied in Thomson 
equilibrium by considering R* as a unit matrix. For the Likins-Pringle equilibrium 
multiply equation (3.73 1) by c;,, and (3.73 ii) by c3, and subtract to get: 
KCNC s5€ 4 C36,,) = 0 (3.74) 
By using the relationships between the elements of R” the expression in the 
parentheses is equal to c;,(Appendix B). So the final form of the equilibrium conditions 
shows that there are two types of Likins-Pringle equilibria conditions: 
Ci). 30S 0 
(i) ¢3,40; c),=0 (3.75) 
For the equation (3.75 i) the spin axis p, is normal to the vertical axis of the orbit 
frame ¢,, but not normal to the forward direction of the orbit frame¢,, so the spin axis 


generates a hyperboloid in inertial space (Figure 3.3). However this condition is not the 
interest of this study. On the other hand the equation (3.75 ii) implies that the spin axis is 
normal to the forward direction but not to the vertical, and the spin axis produces a cone 


in inertial space. The results for conical case can be obtained by solving equations (3.731) 


and (3.73 ii) for c}, #0 and cj, =0. 


(i) 


(ii ) 


From Appendix B: 


Ss AS EIS > Ls 8 
C33 = Cy Cy. — Cpa Coy => 


SS os ° re 
C33 =C11C5) Since c3, =0 


5S Uys .s _ 85,8 S Ls ys . Se J. 
Cry = CC, HCG, “=> e2, =—e7e3,. since cj, = 0 ‘then 


Ss 
s __ © 
C3 = 5 


Cy 


Substituting (3.76 1) and (3.76 11) into (3.73 1) yields: 


Ci SoS Ss Ss 4% S 
k,3(-2 (eye) — Cpe, ) + (I1+k, JU +¥)c3, =0 


11 


k,(-4e2,¢3,)+ +k, + V)c3, =0 


From Appendix B: 
: S 8 ,8 _ 8 ,8 
(1) Cj = C513) — Cy9€3, 


+: S28 8 48 8 
(ii) Cf, = €53€3, — C51C33 


S ps .s = Ss 
=> C3 =—C>,C3, since c,, =0 
SU), 65,6,,.. “siiee > 65, 0 ‘then 


S 
s _ [1 
C3 = 





S 


C3) 


Substituting (3.78 1) and (3.78 11) into the equation (3.73 i1) yields: 


C 





kG 
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2 (—€3,€3,) — Cy C),) + +k, JU +V)ci, = 0 


k,(-4c3,¢),)+ +k, 1 +V)c;, =0 


As aresult, equations (3.77) and (3.79) become: 


(1) 


(11 ) 


C3[—4k,c5, +(1+k, )(1+Vv)] =0 


Ch[-4k,c3, + +k, U+¥)] =0 


(3.76) 


(3.77) 


(3.78) 


(3.79) 


(3.80) 


Since c;, would be zero when both c;, and c;, are zero, these equations can only 
be satisfied when the expression in the parenthesis is zero. The solution for these 
equations is given by: 


(1+k,)+V) 


Ss S An 
c), =0 C5, =cosP®= 
21 22 
4k 


(3.81) 


t 
where ® is the constant rotation angle about the forward axis of the orbit frame. 


The transformation between the stroboscopic frame F’,and the orbital frame F, 


R® =R, R(®) (3.82) 
where ® is the constant tilt of the spin axis about the forward direction of the orbit 
frame and y represents the three small angles which allow small attitude excursions 


about the nominal motion. Therefore the transformation matrix between the stroboscopic 


frame and the orbital frame is: 


1 Vows 92 || 0 0 
Re Saye. oll y, ||0 cos® sin® (3.83) 
Vou 2G: A 0 -—sin® cos® 


By using these relationships, the linearized equations of motion are given by 
Hughes as [17]: 
(i) ¥,-A+k,\(@ sin®)y, -(1+3k,)(@, cos®)y, —4k,(@, sin’) y, =0 
(11 ) ¥,+(@, sin®)y, =0 


(iii) 7,+(1+3k,(@ cos®)y, —3k,@77, =0 (3.84) 
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Because it is impossible to maintain the stability for the spin rate of spinning 


satellite, this research is only concerned with the stability of the spin axis (1.e. y, and 7;). 


After integrating the second equation and substitution in the first equation, new equations 


for the angles y, and y, can be found [17]. 
(i) ¥,-(+3k,)(@, cos®)y, + (1-3k, )(@, sin®)’ y, =0 
(ii) ¥,+(+3k,)(@, cos®)y, —3k,@°7, =0 (3.85) 


The characteristic equations for these two equations are: 


4 2 
[=] +b, (=) +b, =0 (3.86) 

O. Oo. 
b, =14+3k, +9k? —9k (1+kt)sin? ® (3.87) 
b, =-3k,(1—3k,)sin® ® (3.88) 


The necessary and sufficient conditions for infinitesimal stability are: 
b, > 0 
6 >0 


b; —4b, >0 


The regions in which these conditions are satisfied are shown in Figure 3.5. The 
initial design of the TSS will be based on these stable regions from the rigid body 


analysis. 
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RRM RXXX ARAN XXKX XAXK RXXK AXAX RXXX AXKK XXX 
HANK KXAN RXAX AXKX RANK RXXK RXAX AXAX XXXK XXX 
MXRMM RXXM RXAX RMKN RXANM RAR RRAK ARAN RRKK RAK 
MARK RXXX RXXX RXXX RAXX RXXK RXAX ANAK AXKX XXX 

HAN AXA AXRAX AXNK XAXK RXRX RXAX RXAX XAXXK XXX 
MXMM RXXRM RXAN RRKM MXNM RAR ARAN RANMA RRNK RX 
MRNRK RXXK AKAN XXKX RANK RXR AKAN RXAX AXKK XXX 
RANK RXAK AXAX HXKX RAXK RXXK RXAX XXX AXKK XXX 
MMMM RXXM RXNM MRKM MXNM RANK ARAN RRNA RRNK RAK 
MRNK RXXX ARAN AXKX KANN RAK RXAX AXNX AXKK KXX 
MARK RXAK AXAX AXXX KAKA RXXX AXAX AXAX AXXX XXX 
MXMM RXXRM RXXX RRKM MXNN RAR ARAN ARAN RRRK NRX 


MAK RRXX RXAX XX 

HAN AXA AXXX XXX: 

MRM RXXM RARNN HXKK 
MRK RXRX RAXAX XXX XX! 
HAN RXAX AXRAX AXXK XXX 


MRK RXXK RRNK XXKK RXXX XX: 
RAK AXA ARAX HXXK XXXX XXXX 


MRK REX ARAM RXKX RANK RXRK AXRX 
NRK RXAX AXRAK RXXX RXXX XXX XXXX XXX 


HRN RXXK ARXX XXKK XXXX XXX RXXX XXXX XXXX 
RAN AKAX AXAK RANK RXXX RXAX AXXX XXXX RANK KXXX 
RRM RXXK RARNM MXANM MARK RXRKM ARRAN RRNA RANK RRAK 
NRK RXXX ARAN XXKX RXXK RXXK RAXKX KXXX RXXK AXXX 
HAN ARAX ARAX RANK RXXX RXAX AXXX XXXX XXX XXX 
NRX RXRRM RXXM RXXRM RRAN RXNA KANN RARK RRRK RRAX 
MRK RXXX ARAN XXKX RXXX XXX RXXX XXNX RXXK XXX 
RAN RXAX ARAX AXKK RANA KXXX RXAX AXNX XXX XXX: 
XXX RXRMM RXXM RXAM RRAN RAXNA RRANM RANK ARAN ARKX 
MRK RAN ARAN XXKX RXXK XXX RAXKX KXXX RXXK RXXX 
RAK RXAX ARAK AXAKK RXXX AXAXX RAXXX XXXK RXXK XXX 
XXX RXRMRM RXXM RXRM RRNN RNNA MXANK RXRK ARAN ARKX 
MRK REN ARAN KXNX RANK RXXX RXXX RXKX RXXK RXXX 
MAK ARAN ARXX AXKX RXXX ARAN AXXX XXXX RXXK RXXX 
NRX RXRMRM RXXM RXAN RRM RANMA RARM RXRK RRAK RRAK 





Figure 3.5 Stability Diagram for the Likins-Pringle Equilibria [17] 
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3.2 Design of a Tethered Satellite System 


The tethered satellite system which is the subject of this thesis consists of three 
satellites at the corners of an equilateral triangle with the same mass, and two end bodies 
at the opposite sides of this triangular plane. Satellites and end bodies masses are 


considered as point masses. The tethered satellite system is shown in Figure 3.6 


1. end body 





2. end body 


Figure 3.6 Design Of The tethered Satellite System 


In Figure 3.6 mass of the satellites are given by m, and the mass of the end bodies are 
given by the termm,. D gives the distance between the satellites in equilateral plane, and 


L is the distance between the two end bodies. Because of the axisymmetry of the 
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formation, the center of mass of the system is on the triangular plane and at the center of 
the this triangular plane. Position of the satellites and the end bodies with respect to the 
center of mass of the system are: 

(a) 7, =[0 -L/2 0] 

(b) 7 =[0 1/2 0 

(c) F=Lp/2 0 DiQv3)| 

(d) F=l0 0 -p/V3| 

(e) F=[D/2 0 DIQ2v3)| (3.89) 
and the total inertias of the system is: 


Palit tel atl (3.90) 
where the terms /,, and /,, are the inertias of the end bodies, and the terms /,,, /,,, and 


I,, indicate the inertias of the satellites 











lm. D*m 
—— + . 0 0 
f= 0 D’m, 0 =|0 J, 0 (3.91) 
2 2 
0 0 mb 2 mM, 0 0 7, 
2 2 | 








where the axial and transverse inertias are J, and J, respectively. The variable k, is 


defined as: 





k, =—! (3.92) 
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Substituting in the inertias for the TSS yields: 


2: 2 
Pi m,L ge m, 
2 2 


mL’? d’m, 
+ 
2 2 








k,== (3.93) 











To analyze the parameter k,, it is possible to define non-dimensional variables: 


d= and m=— (3.94) 


m 


e 





D m 
L 


Therefore equation (3.93), can be expressed as: 





(3.95) 
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Figure 3.7 k, Values for Mass and Size Ratios 
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The Likins-Pringle relative equilibria states that there is a constant rotation angle 
about the forward direction (¢,) of the orbit. The spin axis of the satellite generates a 
cone in inertial space [17], and the rotation angle about the forward direction gives the 
half cone angle. The equation for the half cone angle is: 


(1+k,)1+V) 
4k 


Cy, =cosP= (3.96) 


t 


fat and P= (3.97) 





where: k= 


Reference frames for this orientation and the half cone angle is shown in Figure 3.8. 


Note that the spin axis of the formation is $,, so a half cone angle of 90° implies the 


formation is aligned with the local vertical. 





Figure 3.8 Reference Frames 
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To analyze the effects of the spin rate and the distance parameters of the TSS 
design on the equilibrium position of the formation, equation (3.96) can be used. Figure 


3.8 shows the half-cone angle for different L/D and spin ratios with constant mass ratio of 


m,/m,=2. 


Half Cone Angle 





4 6 8 10 12 14 16 18 20 
L/D 


Figure 3.9 Half-Cone Angle for Different Spin Rates 


In Figure 3.9 the half cone angle increases with increasing positive spin rate V or 
decreasing size ratio L/D. To obtain the equilibrium position of the tethered formation 
at nadir direction the half-cone angle should be 90° deg, which is achievable with the spin 


rate: 


(3.97) 


_— 
II 
| 

— 
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According to equation (3.56), this spin rate implies that the satellite doesn’t spin 
about its symmetry axis with respect to the stroboscopic frame. However a tethered 
satellite system, which consists of distributed satellites must rotate about the symmetry 
axis to rigidize the system and to obtain the stability. The spin rate about the symmetry 
axis of the tethered formation with respect to the stroboscopic frame is: 


py AVEO, 





s (3.98) 
Qo. 
Substituting V into the equation (3.98) yields: 
V,2V+1 (3.99) 


By using this relationship it can be seen that there is symmetry in the Figure 3.8 
for the: 


yp. =0 (3.100) 


This result implies that positive spin rates about the symmetry axis b, in the body 
frame gives the half-cone angle greater than 90° degree and negative spin rates about the 
symmetry axis b, in the body frame gives the half-cone angle less than 90° degree. 

The other way to get a half-cone angle closer to 90 degree is to increase the size 
ratio of formation, which is denoted by L/D in Figure 3.9. Increasing the distance 
between the two end bodies increases the ratio L/D and causes to the half-cone angle to 
approach 90°, which implies the spin axis is on the nadir direction. 

To analyze the effects of all parameters of the TSS design, (spin rate, mass ratio 
and the size ratio) on the equilibrium position of the formation on the orbit, non— 


dimensional variables in equation 3.94 and equation 3.95 can be used. By using the non- 


3IR 


dimensional variables k, and the spin rate V, the rotation angle about the forward orbit 


direction is shown in Figure 3.10. 
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Figure 3.10 Half-Cone Angle for Conical Stability 


For the Figure 3.10 there is symmetry for the half cone angle about the spin rate: 
v=-1 (3.101) 
As with the analysis before, this implies that the satellite does not spin with respect to 
the stroboscopic frame. To design a tethered formation in the stability region, which has 
the spin axis closer to the nadir direction, the spin rate should be as small as possible and 


the variable k, must be negative and close to —1. 
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IV. Analysis and Results 


To analyze the dynamics of the tethered satellite system, equations of motion for 
each satellite are developed and solved numerically. In this chapter, forces on the 
satellites and the end bodies will be determined. In this study, perturbation forces like 
earth oblateness, gravity effect of sun or the gravity effect of moon are not considered. 

To keep the tethered satellites at the desired positions with respect to the each 
other, the system must rotate at some minimum rate. Parameters that affect the minimum 
spin rate will be explained below. 

Numerical solutions for the equations of motion will be explained using the 
Likins-Pringle equilibria. The stability region will be presented where the tethered 
satellite system is stable (by the definition of the stability used for this study). Specific 
examples in this region will be explained in more detail. Stability of the systems for these 
examples will be discussed and the parameters which have effects on the stability of the 


systems will be explained. 


4.1. Equations of Motion 


The equation of motion for a satellite in a inverse square gravity field is [19]: 


mr =———7 (4.1) 
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where is the earth gravitational parameters, m is the mass of the satellite and 7 denotes 


the position of the satellite from the center of the earth. However, due to the perturbing 
forces transmitted to the satellites through the tethers, the motion of each satellite in a 
tethered system is non-Keplerian. For a simple tethered satellite system, where the tether 
is considered as a massless and extensible string, forces on the particles can be shown in 


Figure 4.1 [18]. 





Figure 4.1 Three Satellites With Tethers 


Three satellites, which are connected to each other with tethers are shown in 
Figure 4.1. Their distances from the center of the earth are presented with 7, vectors. 
There are two types of forces on each satellite, which are the gravitational force and the 
tension force. Gravitational force depends on the distance of satellite from the center of 
the earth, and tension force depends on the distance between satellites and the tether 


length. Equations of motion for each satellite are: 


4.2 


= HM, — No 
ny heh a (4.2) 
Al i. 
= Hm, — Np 9 
ze aa ie | os (4.3) 
7, | FoI 
= Hm, — T9 
WT tam) ee (4.4) 
| | 


where m,is the mass of the satellite and, 7; is the tension force between bodies i and j 


and: 


ho=%-7 and 7,=7,-7, (4.5) 


r,. 
The — terms show the direction of the tension force. 


Vj 








Magnitude of the tension force can be given by [7]: 


Re E(j-FI-1,) #1, ¢ 
. ely 





Fade 


(4.6) 





7, 
where F is the tether extension stiffness and /, is the unstretched tether length. It is 


obvious that the tension force occurs only when the distance between the satellites is 
greater than its initial value and there is no compressive force. 

With these basic equations, equations of motion for the Tethered Satellite System 
model, which consists of three satellites at the corners of an equilateral triangle and two 
end bodies at opposite side of this plane, connected with massless and elastic tether to 
each other can be presented. Figure 4.2 shows the orientation of the tethered satellite 


system on orbit. 
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Figure 4.2 Orientation of TSS on Orbit 
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(4.8) 


(4.9) 


(4.10) 


(4.11) 


4.2 Minimum Spin Rate For Stability 


For a rigid satellite formation system, rotation of the system can be used for 
attitude control. On the other hand, the tethered satellite system with flexible tethers must 
rotate to keep the satellites in the formation at their nominal position with respect to each 
other. The TSS considered for this study consists of three satellites and two end bodies. 
The three satellites are placed on a triangular plane. With some simple analysis it is 
possible to determine the minimum spin rate to keep the formation from collapsing. 

The aim here is to keep the triangular plane of three satellites perpendicular to the 
nadir direction as much as possible. It is not possible to keep the satellites on a planar 
shape with a definite distance to each other without spinning. Minimum spin rate for the 
tethered satellite system can be found by assuming the spin axis of the tethered system is 
aligned with the nadir direction. 

For the TSS model, the angular velocity of the system normal to the orbital plane 


can be calculated from the orbital mean motion: 


QO = 


Cc 


yu 
S (4.12) 


The centrifugal force is: 
F.=mra: (4.13) 
Thus for the TSS model the net force on the end bodies in the formation is: 


sr add a (4.14) 
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C, 1.end body 


Figure 4.3 External Forces on the first end body 


External forces on the first end body are presented in Figure 4.3, and by using the 


distance of the end body from the center of the earth: 


Neb = (Fon -L/2) é; 


(4.15) 


where r,, denotes the distance of the center of mass of the system from the center of the 


earth. Gravitational force on the first end body is: 


_ em, - 
3 Topi 
Vepy 


gl 








where m, is the mass of the end body. 
Substitute the equation (4.15) into the equation (4.16): 


nl Hm, a 
FS eG De, 
(r,, -L/2) 





The centrifugal force on the particle is: 


= fas 2 
as =MN Fp) @ 


Cc Cc 
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(4.16) 


(4.17) 


(4.18) 


Substituting the equations (4.12) and (4.15) into the equation (4.18) the centrifugal force 


on the end body is: 


F =m £12) —<, (4.20) 


Thus the net external force is: 





Fen ie ay —L/2)eé, ay Coes ye ead (4.21) 
(Pa, — LE /2)| an 
ie 
_ ics 
T, 


183 


Figure 4.4 Net Forces on The End Body 
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Figure 4.5 Tension Force On A Tether 
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The total force on the first end body and the tension force on the tether are 


presented in Figure 4.4 and Figure 4.5. For these figures 7,,denotes the tension force 


1si 
between the first end body and the satellite. Z is the distance of the end body from the 
center of the mass of the system and Dis the distance of the satellites from the center of 
mass. The angle between the symmetry axis and the tether is: 


D/N3 
D 


@=arctan 





(4.22) 


Because of the axisymmetry of the TSS model tension force on the tethers 


attached to the end body are approximately equal to each other in magnitude, so: 


| eT (4.23) 


net| 2 1si 


Thus the tension force on the tether is: 





=e ae 
T,, =. —_ 4.24 
3° cos Con 
where 7,,, denotes the tension force on the tether between the first end body and the i th 
satellite. 


By using the same method for the second end body, which is at the opposite side of 
the triangular plane, the total force can be calculated. Position of the second end body is: 
hoe (FPL 2), (4.25) 
Gravitational force on the first end body is: 


Mm, _ 
alt aan (4.26) 


g2 








Pop? 


Substituting the equation (4.25) into the equation (4.26): 
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F 5 +L/2)é, (4.27) 
(r,, +L 12) 





The centrifugal force on the particle is: 


= = 2, 
FP spa =F 59 @ 


Cc Cc 


(4.28) 
Substituting equations (4.12) and (4.25) into equation (4.28) the centrifugal force 

on the end body is: 
F 


cle 


=m,(r,,+L/2) © é, (4.29) 


cm 3 
cm 


Thus the net external force is: 


a ae OME, HON MES OE ie) eer (4.30) 
(Pon +L /2)| r, 





the geometry is similar to the first end body, so the tension force on the tethers 
between the second end body and the satellites are: 
es oe | 
= 4.31 
*“ 3 cos oe 





The forces on the satellites, which are placed on the equilateral triangular plane, 
are the tether tension forces, the centrifugal force due to the spinning of the formation, 


and the gravitational force. 
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Figure 4.6 Forces On The Satellites 
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Where 7. 


,, and 7,,,denote the tension forces, F’,,, denotes the component of the 


gzsi 


gravitational force in z direction and F.,, denotes the centrifugal force on the satellites. 


Here F.,, is the centrifugal force due to the spinning of the satellite. The component of 


the angular velocity due to the rotation of the formation around the earth is not included 
in this analysis since it will vary depending upon the location of the satellites. 


The total force on the satellite is: 





Fa? eet dla (4.32) 


tot gzsi 1si 2si 
The centrifugal force on the satellites due to the spinning around the axis 


perpendicular to the triangular plane is: 
Fj.=m,—v" (4.33) 


where the term Vv, denotes the spin rate about the axis perpendicular to the plane 
with respect to the stroboscopic frame, and m, is the mass of the satellites. 
V.=V+Q., (4.34) 


To get the equilibrium spin rate, set F,,, in equation (4.32) equal to zero and solve 


tot 


for V,: 


eS c= a= 
Ve = a thy +T,,,) (4.35) 


Ss 


This is the minimum spin rate to keep the tethers in tension. Figure 4.7 shows the 


values of the v, fora range of L/D and m,/m,. 
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Figure 4.7 Minimum Spin Rate For Stability 


Figure 4.7 shows that the minimum spin rate increases with the mass ratio, which 
is the ratio of the mass of the end body to the mass of the satellite for this figure. Also, 


the minimum spin rate doesn’t depend on the distance between the end bodies. For the 
mass ratio 0.5 the minimum spin rate is approximately the J2 times greater than the 


nominal motion @.. 
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4.3 Simulation Results 


The possibility of using Likins-Pringle relative equilibria for the design of the 
tethered satellite formation is demonstrated using numerical simulation. The equations of 
motion have been solved numerically because of their complexity. 

The example provided here is a tethered satellite formation in a circular orbit with 
a radius of 8000 km. The formation is shown in Figure 4.8. Tether extension stiffness is 


assumed as 10° and constant for all simulations. 
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Figure 4.8 Reference Tethered Formation in the Body Frame 
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The formations parameters are: 


m =m 


satl 


M erod\ ca MN ehody2 = m, =100kg 


sat2 ~ """sat3 


=m,=200kg (Mass of the satellites) 


(Mass of the end bodies) 


L=80km (distance between two end bodies) 
D=10km (distance between satellites) 
o,= ae (nominal motion of the tethered formation) 
Vom 
V,=-2*o, (spin rate with respect to the stroboscopic frame) 


Transverse and axial inertias of the formation are: 


Ef “Pp 
me 4 = Ms _ 330,000 





I, =D°*m, =20,000 
Substituting the spin ratev, into equation (4.34): 
V=V,-@,=-3* @, 
Then the spin ratio is: 


vat =-3 


V 
QO. 
The inertia ratio k, is: 


k == =— 0.9394 


t 





t 


(4.36) 


(4.37) 


(4.38) 


(4.39) 


(4.40) 


For the equilibrium position of the tethered formation the half cone angle is: 


(1+k,)+V) 


= ateood seis deg 


t 


4.13 


(4.41) 


The angle between the spin axis and the nadir direction is: 
90 —& = 1.8486 
Figure 4.9 shows the orientation of the reference tethered formation in the orbital frame 


with an angle of 1.8486 degree between the spin axis and the nadir direction. 





Figure 4.9 Tethered formation on the Orbital frame 


With this orientation the satellites follow a circular path because of the formation 
spin, and the projection of the satellites on the Earth surface will be approximately a 
circle. Figure 4.10 shows the satellites and end bodies during 60 orbits in the orbital 
frame and Figure 4.11 shows the formation on the s)-s3 plane of the stroboscopic frame. 
Tethers between the satellites and the end bodies are not shown to simplify the figure. 
Simulation results of this configuration show that the tethered formation is stable up to 60 


orbits. After 60 orbits the satellites in the formation depart from their initial position. 
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Figure 4.10 Formation During 60 Orbits in Orbital Frame 
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Figure 4.11 Positions of the Satellites and the End Bodies 


in the Stroboscopic Frame 
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4.4. Analysis of the Design Parameters 
In this section the parameters that affect the stability of the tethered formation will 
be investigated. Stability of the tethered formation is affected by: 
Spin rate of the formation 
Mass ratio 


Size of the formation 


4.4.1. Spin Rate 
Spin rate of the tethered formation changes the rotation angle, ®, about the ¢, 
axis. Increasing spin rate would decrease the rotation angle and the tethered formation 
will be less stable. Effects of the spin rate on the stability of the formation are 
investigated here for a different spin rate. The same mass and the size ratio of the 
reference configuration have been used. 


For the spin rate v, = -3* @., half cone angle will be ® = 87.2265. This spin rate 


is higher than the spin rate, which is presented as a reference in the previous section. This 
configuration of the tethered formation is not stable for long period of time, and after the 
second orbit tethered formation begins to depart from the desired configuration. Figure 
4.12 shows the positions of the satellites and the end bodies after the second orbit in 
orbital frame. Figure 4.13 shows the positions of the satellites and the end bodies in the 
stroboscopic frame for the same orbits. 

This simulation demonstrates that the stability of the tethered formation decreases 


with the increasing spin rate with respect to the stroboscopic frame. 
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Figure 4.12 Positions of the Satellites and the End Bodies in 


The Orbital Frame for Higher Spin Rate 
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Figure 4.13 Positions of the Satellites and the End Bodies in 


The Stroboscopic Frame for Higher Spin Rate 
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For this configuration, satellites and the end bodies begin to depart from the initial 
position in the second orbit and for long period this deviation will be worse. This 
behavior of the tethered formation can be explained by the reduced half cone angle. The 
formation spin axis is farther from the nadir direction, so gravity gradient torques 
destabilize the system. Figure 4.14 shows the configuration of the tethered formation 


after 10 orbits. 
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Figure 4.14 Tethered Formation with Higher Spin Rate 
After 10 Orbits in Body Frame 
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4.4.2 Mass Ratio 
The mass ratio of the tethered formation changes the minimum spin rate of the 


configurations and the inertia ratio, k,. Different values for k, cause a different half cone 


angle and orientation on the orbit. Effects of the mass ratio on the stability of the 
formation are investigated here for a different mass ratio. The same spin rate and the size 
ratio of the reference configuration have been used. 


For the mass ratio m,/m, =3, half cone angle will be ® = 87.1810. This mass 
ratio of m,/m, is higher than the mass ratio, which is presented as a reference in the 


previous section. The tethered configuration is not stable for long period of time. After 
the fifth orbit, the tethered formation begins to depart from the desired configuration. 
Figure 4.15 shows the positions of the satellites and the end bodies during eighth orbit in 
orbital frame. Figure 4.16 shows the positions of the satellites and the end bodies in 
stroboscopic frame. 

This simulation demonstrates that the stability of the tethered formation 


decreases with the increasing mass ratio between the satellites and the end bodies. 
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Figure 4.15 Positions of the Satellites and the End Bodies in the 


Orbital Frame for Higher Mass Ratio 
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Figure 4.16 Positions of the Satellites and the End Bodies in 


The Stroboscopic Frame for Higher Mass Ratio 
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Figure 4.17 shows the tethered satellite formation configuration after 10 orbits. 





Figure 4.17 Tethered Formation with Higher Mass Ratio 
After 20 Orbits in Body Frame 


4.4.3. Size of the Formation 


The size ratio of the tethered formation changes the value of the inertia ratio, k,. 
Different values for k, will cause a different half cone angle and orientation in the orbit. 


Effects of the size ratio on the stability of the formation are investigated here for a 
different size ratio. The same spin rate and the mass ratio of the reference configuration 


have been used. 
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For the size ratio L/D=9, half cone angle will be ® = 88.5493. This size ratio of 
L/D is higher than the size ratio, which is presented as a reference in the previous section. 
Figure 4.18 shows the positions of the satellites and the end bodies for three orbits 
between the 65 and the 68 orbits in the orbital frame. Figure 4.19 shows the positions of 
the satellites and the end bodies in the stroboscopic frame for the same orbits. After 68 
orbits the shape of the formation begins to deform. This simulation demonstrates that the 
half cone angle increases with increasing size ratio and the angle between the spin axis 
and the nadir direction decreases with higher size ratio. Therefore increasing the size ratio 


increases the stability of the tethered formation. 
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Figure 4.18 Positions of the Satellites and the End Bodies in the 
Orbital Frame for Higher Size Ratio 


4.22 







Wo 





£S 1.end body 
O 2.end body 
© 1.satellite 
© 2.satellite 







ty See eee Cee 









s3 
o 
: 
: 
' 
' 
: 
: 
eee ee ee 
: 
: 
: 
; 
: 
; 
: 
: 
H 
: 
: 
' 
: 
; 
; 
B 
' 
: 
; 
: 
; 
' 
; 





eee eee eee 








Ee ob connate 
“Ra eee 






Figure 4.19 Positions of the Satellites and the End Bodies in the 
Stroboscopic Frame for Higher Size Ratio 


Figure 4.20 shows the tethered satellite formation configuration after 80 orbits. 





Figure 4.20 Tethered Formation with Higher Size Ratio 
After 80 Orbits in Body Frame 
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4.5 Stable Configurations 


Simulation results show that the stability period of the tethered satellite formation 
depends on the configuration of the system. The half cone angle, which is the result of the 
design parameters, determines the stability of the tethered formation. Because of the 
gravity gradient force, configurations, which have the spin axis closer to the nadir 
direction, are more stable than the other configurations. The tethered formation which has 
the spin axis aligned with the nadir direction would be stable for very long time. 
However, to maintain such configuration is not possible for the non-rigid tethered 
satellite formation because of the minimum spin rate constraint. 

Analysis shows that the stability decreases dramatically for the tethered formation 
configurations which have a half cone angle smaller than 88 degrees. Because of this a 
simulation has been set up to find the stable configurations with a minimum 88 degrees 
half cone angle; that is, the maximum deviation of the spin axis from the nadir direction 
is allowed to be 2 degrees. Figure 4.21 shows the stable configurations of tethered 
satellite formation with the half cone angle between 88°and 90° degrees for the spin rate 


with respect to the stroboscopic frame on the order of the orbital rate, which is: 





. V+a, 
VS 
o 


Cc 


(4.42) 


and the size ratio of L/D in the range between 6 and 15. This plot supports the 
conclusions of the previous sections. The stability increases with larger L/D and lower 


spin rates (up to the minimum spin rate). 
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Figure 4.21 Stable Configurations 
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V. Conclusions and Recommendations 


Likins-Pringle relative equilibrium has been analyzed for the possibility of using 
this equilibrium to find the stable configurations of the tethered formation about near- 
nadir directions. Since the application of the tethered formation is Earth surveillance and 
sensing, the system was designed so the sub satellites are widely distributed about the 
nadir directions. Two end bodies, which are placed at the opposite side of the triangular 
plane help for the system stability by using the gravity gradient force. 

Stability of the tethered satellite system depends on several design parameters. 
Obtaining a stable tethered satellite formation with a spin axis perfectly aligned with the 
nadir direction is not possible. However it is possible to obtain a stable formation system, 
which has an angle between the spin axis and the nadir direction from very small values 
to about 2 degrees. 

For this study the tethers between the satellites and the end bodies are considered 
as massless, extensible tethers. For the future study, the tethers should be modeled with 
one of the discritization methods to more accurately model the flexible nature of the 
tethers. 

For this study the effects of the air drag, moon’s gravitational force or the 
gravitational force of the sun are not considered. For the future study about the tethered 
satellite formation these perturbation forces should be taken into account and the effects 


of these forces on the stability of the tethered formation should be demonstrated. 
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Appendix A 


Linearized Solution of the Equations of Motion 


Equations of motion in the linearized form: 
(i) & +(-k,)va, +k,v°Q, +3@7k,(a,c+a,s)c =0 
(ii ) a, =0 
(iii) & —(—-k, va, +k,v’a, +3@7k, (a,c + a,s)s =0 


where, 


s=sin(v+@.)t c=cos(V + @. )t O, = 3 





(A.1) 


and 


(A.2) 


The linearized equations of motion include time variable coefficient because of 


the nominal motion of the satellite and the reference spin rate. It is possible to simplify 


these equations by transforming them to a new set of coordinates that will eliminate the 


periodic coefficients. The new parameters are defined as [17]: 
(i) 4¥,(@¢)=a@,(t)cos(v +a, )t+a@,(t)sin(v + a. )t 
(li) 7% (t)=@,(t) 


(ii) ¥,()=-a@,(t)sin(v+ @,)t+ @,(t)cos(v + @, )t 


(A.3) 


Equation (A.1 i) implies that the satellite is unstable about the spin axis. Thus for 


the stability about the other two axis the equations of motion can be written in matrix 


; le 0 eae bl 
+ + 
0 1]|/a@,| |-d-k,v 0 a, 


3a@7k,c? +k,v? 307k, sc a,|_|0 
3a@- k,sc 3a@2-k s’ +kv? |La,} [0 


And the new set of coordinates: 


Pee i 


Inverting and taking derivatives of this equation gives @,,@,,@,, @,,@,,@, in 


form as: 


(A.4) 


terms of the new set of coordinates 7,,7;,7,,73,);.7;: 


eerie ba ie 


By taking the derivative of the equation (A.6): 
lvl cles fn] = 
a, s chy, 6: S596 


p=(V+@. ) (A.8) 





where, 


And the second derivative of this equation is: 


alls clade alee cll Es “le 
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By adding the similar terms: 
Sk Bet BEE) « 
a, SC LY; c —SILYs —S —CTVWY; 


er a) (a | |al. 
Substituting the equations for, ‘|| il ie into the equations of motion 
a 





Pedals walle ls la) 
S  ¢€ ILY; Cc —SILYs —S CTL; 
0 rid C ple . | 0 ula 
_ |+p + 
—(1-k,)v 0 Is € 1% c —s ||-(l-k,)v 0 Va 
ee +kyv iia : é a )-14) (A.10) 
307k ,sc 30 ;k,s  +kv |\is c || 7; 0 
Oe ede seared 
8 36) Ge 8 cts? 0 1 0 
8 || fer gO" al (A.12) 
c -s| ls c]l1 0 
Substituting the equation (A.12) into the equation (A.10): 
FAR Eee EMay P UE) 
Se y, 1 O% y,| |-d-4,)v 0 Vs, 
ok | 0 a ene 0 (Hl ae 
1 0 ||-d-k,v 0 Y; 0 kv? || 7; 0 








Note that: 


Thus: 
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Simplifying the equation (A.13): 


|: mal 0 areca. 
SC V; 2p-U=«,) 0 V3 


—p’ + p(l-k, v+3a@2k, +kv 0 Va 0 
p + p-k,) E Wits (A.14) 
0 —p +p(l-k,)vt+k,v" | 7; 0 


Since: 





Cc. = 
ts] |-120 (A.15) 
S Cc 


The term in the parenthesis in the equation (A.14) must be equal to zero to satisfy 


the equation: 


aul 0 awa Mar 
V; 2p AEK,) 0 1 


eed _ 2 0 
p + pd-k,)v+307k, +k,v 0 NM \(_ (A.15) 
0 —p° +p(l-k,vt+kv" || 7% 0 


Substituting the equation (A.8) into the equation (A.15): 
p=(V+@,) 
Elements of the equation become: 


—2p+(1-k,)v=-2(v+a.)+(1—-k,)v 


(A.16) 
=-(V+2@,+k,V) 
2p-(1-k, v=2(v+a,)—-(1-k, )v 
p—(-k, v=2( Uk) (A.17) 
=Vv+20.+k,v 
—p? + pl-k, v+3a@2k, +k v=-(v+,) +(V+a,)1-k, v+30@2k, +kv 
=-VO, —@. —k,vo, +3k,@2 A.18) 


=k, (30; —v@,)-@,(V+@,) 
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—p’+pl-k, vw+kv’? =-(Vt+a,) +(v+a,)\1-k, vtkv 
=-VO, —@; —k,vo, (A.19) 
=-(k,vo,+0(V+@.,)) 


Substitute the equations (A.16) through (A.19) into the equation (A.15), equation 


of motion in the matrix form becomes: 


y, - 0 —(V+20,+k,v) || " 
y; (v+20, +k,v) 0 y; 


(A.20) 
k, (3; —vo,)-@,(V+@,) 0 y,)_ [0 
0 % (k,Vvo, oF QO. (v ae Q. )) 3 0 
Finally the equations of motion can be written as: 
(i) ¥-V+20, +kV)y, + (k, 30; —VO,)-@,(V+@,))7, =0 
(11 ) 7, =0 
(il) ¥,+(V+20,+k,v)y, -—(k,vo, +@,(V+@,))7; =0 (A.21) 
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Appendix B 


Rotation Matrix 


All vectors represented by their components with respect to the reference frame. 
This representation is relative to the selected reference frame, and it is possible to write 
the same vector with respect to another reference frame. To write the vector in another 
reference frame, transformation is used between these reference frames. Rotation matrix 
gives the transformation of any reference frame to another. Elements of the rotation 
matrix show the relationship between the axes of two reference frames, and they are the 
cosine of the angles between each axes of the reference frames. Because of this, rotation 


matrix 1s also called as the direction cosine matrix [18]. 


a, ,b 





Figure B.1 Rotation of the Reference Frame 
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f° 2750, +1 5a, +30, (B.1) 
The term 7“ is the representation of r vector in F’, frame. 

aed Fear, 8 sae tt Ae (B.2) 
The term 7” is the representation of 7 vector in F, frame. Transformation between the 


F' frame and fF, frame can be shown by using the rotation matrix: 


F,=R™F, where R”’is the rotation matrix R™ =|c,, Cy ©; (B.3) 


> 


b, Cy Cyn O43 ay 
by |=] Cy, Co. Ca3 | | (B.4) 
b, C31 C30 C33 a, 


From the equation above, the relation between these two frames F’, and F, can be 
written for each b, ; b, and b, [17]: 
(i) b, = C114, + C24, + C34, 
(ii) b, = Cy), tly Gd, +Cy,4, 
(iii) b, =c,4, +64, +044, (B.5) 


By using the cross product rule equations for the relation between the axis can be 


written: 


An 


Gd, X ay = a; DO Gy Gq CGS a5 (B.6) 


An nA 


b, x b, = b, Dy & b, =b, b, xb, = b, (B.7) 


= 


From the equations (B.5): 
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A 


, 


A 


b, = 


A 


be 


Dy = (Cy), +. Cy9@, +C 533) X (C314, +034, +0534) 


= (Cy C93 — CygCyy Oy + (CoxC5y — Cp O33 )Og + (Cail ap — Cop As 


b,=b, xb 


An 


Dy = (€5)@y ACy94y-+ C5345) X (CG, FCpGy + C445) 
(C5913 — Cy3 C19 )Ay + (C3301) — Cq4Cy3 By + (C51 Cyy — Cy9€y) AS 
b,=5, xb, 


by = (€),4, + C94, +.C,343) X (Cy14, +Cy)G, +5343) 


(C1p€ 95 — C1gC ap JG, + (yg ny — C403 Ey + (Cylon — Cyn Aa 


(B.8) 


(B.9) 


(B.10) 


With the equations above, nine relationships between the elements of rotation 


matrix can be defined [17]: 


Cy = Cy9€33 — 43039 


Cyy = C3031 — C1033 


Coy = €39C13 — €33C 19 


Coy = C3301) — C3133 


C13 = €y1C39 — Cy C31 (B.11) C43 = C31 C yy — C39] 


C31 = Cy90n3 — C130 a9 
C39 = Cp3C91 — C1 1093 


C33 = Cy1Co. — C29] (B.13) 
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(B.12) 
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